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Abstract 

We apply the 3D reduction method we recently proposed for the N-particle 
Bethe-Salpeter equation to the 4-particle case. We find that the writing of 
the A'' > 4 Bethe-Salpeter equation is not a straightforward task, owing to 
the presence of mutually unconnected interactions, which could lead to an 
overcounting of some diagrams in the resulting full propagator. We over- 
come this difficulty in the A'^ = 4 case by including three counterterms in 
the Bethe-Salpeter kernel. The application of our 3D reduction method 
to the resulting Bethe-Salpeter equation suggests us a modified 3D re- 
duction method, which gives directly the 3D potential, without the need 
of writing the Bethe-Salpeter kernel explicitly. The modified reduction 
method is usable for all A'^. 

PACS ll.lO.Qr Relativistic wave equations. 

PACS 11. 10. St Bound and unstable states; Bethe-Salpeter equations. 
PACS 12.20.Ds Specific calculations and limits of quantum electrodynamics. 

Keywords: Bethe-Salpeter equations. Salpeter's equation. Breit's equation. 
Relativistic bound states. Relativistic wave equations. 
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1 Introduction. 

The Bethe-Salpeter equation || is the usual tool for computing relativistic 
bound states [1-20]. The principal difficulty of this equation comes from the 
presence of N-1 (for N particles) unphysical degrees of freedom: the relative 
time-energy degrees of freedom. In a recent work we built a 3D reduction 
of the two-fermion Bethe-Salpeter equation around an unspecified instantaneous 
approximation of the Bethe-Salpeter kernel, the difference with the exact kernel 
being corrected in a series of higher-order contributions to the final 3D po- 
tential. This potential being not explicitly hermitian, we performed a second 
symmetrizing series expansion at the 3D level. The result of the combination 
of both series gave us an hermitian potential which turned out to be indepen- 
dent of the starting instantaneous approximation, and was in fact a compact 
expression of the potential that Phillips and Wallace jl^ compute order by or- 
der. We found that this 3D potential was also obtainable directly by starting 
with an approximation of the free propagator, based on integrals in the relative 
energies instead of the more usual (5— constraint (integrating propagator-based 
reduction). In this form, the method was easily generalizable to a system of N 
particles, consisting of any combination of bosons and fermions. Special cases 
of N=2 or 3-fermion systems were examined. Taking the retarded part of the 
full propagator at equal times, following Logunov and Tavkhelidze (for two 
particles) and Kvinikhidze and Stoyanov ||2l| (for three particles) also leads to 
the same 3D equations. 

For more detailed explanations about our integrating propagator-based re- 
duction, about related problems like Lorentz covariance, cluster separability or 
continuum dissolution, or for comparisons with other approaches, we refer to 
our preceding work [ pO[ . 

In the present work we test our reduction method in the four-body problem. 
The first step is of course the writing of the starting Bethe-Salpeter equation 
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itself. It is an homogeneous equation for a Bethe-Salpeter amplitude, which is 
derived by identifying the residues of the bound state poles in the inhomoge- 
neous Bethe-Salpeter equation for the full propagator G. The interactions are 
introduced via the Bethe-Salpeter kernel K, which is such that the inhomo- 
geneous Bethe-Salpeter equation reproduces, by iterations, the full propagator 
G as deduced by the Feynman's graphs method. The writing of the Bethe- 
Salpeter kernel K is straightforward in the two and three-body problems. In 
the four-body problem, however, we meet a new difficulty due to the mutual 
unconnectedness of the three pairs of two-body irreducible kernels (12) (34), 
(13)(24) and (14)(23): adding simply the six two-body irreducible kernels (plus 
the four irreducible three-body kernels and the irreducible four-body kernel) 
would lead to an overcounting of some graphs in the expansion of G. We would 
get for example the sequences of kernels (12) (34) and (34) (12), which are in fact 
a unique graph or part of graph (unlike sequences like (12) (23) and (23) (12) 
which are different). It is however not very difficult to correct this overcount- 
ing by including three counterterms in K. Computing then the first terms of 
the 3D potential by our 3D reduction method, we find that, at the often used 
approximation in which we keep only the two-body kernels, the (relative and 
total) energy-independent part of these kernels and the positive-energy part of 
the propagators, the final 3D potential is simply the sum of the six two-particle 
potentials. 

It is clear that the difficulty of the overcounting of the sequences of un- 
connected kernels will not be so easily overcomed for > 5. This suggested 
us to modify our 3D reduction method in order to avoid the explicit writing 
of the Bethe-Salpeter kernel K. Our 3D potential is now written directly in 
a straightforward way for any number of particles, with simply the additional 
prescription of removing the duplicating diagrams which appear when A'^>4. 

In section 2, we recall our 3D reduction method of ||2^ for the N-body Bethe- 
Salpeter equation. In section 3 we apply it to the four-body problem. In section 
4 we suggest a general 3D reduction method starting directly from the expansion 
of the full propagator G, as given by Feynman graphs, without the need of 
writing explicitly the kernel K of the Bethe-Salpeter equation. Section 5 is 
devoted to conclusions. 

2 Integrating propagator-based reduction for the 
N-body Bethe-Salpeter equation. 

The inhomogeneous and homogeneous Bethe-Salpeter equations can be written 

G = G" + G"A'G (1) 
$ = G"K^ (2) 
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where G is the full propagator and $ the Bethe-Salpeter amplitude. The free 
propagator G'^ for a system of / = 0, l...iV fermions and b = N-f bosons is 

U — ... Lrj Ltj_|_]^ ... LtjVj K'^) 

the propagators of the fermion i and the boson j being respectively 

T— — (4) 

Pio - hi+ lehi 

G° = ^ ^ J_ V (5) 

where pi is the 4-momentum of particle i , and 

= di .pi + PiTrii, Ei ^ {pf + mj)^ , (jj = ±1. (6) 

We do not specify the reference frame in which we write noncovariant quantities 
like hi or Ei. Our 3D reduction will in fact be frame-dependent. Practically, 
we shall choose the global rest frame. The kernel K will be chosen in such a 
way that equation (|l|) gives by iterations the usual expansion of G in terms of 
Feynman graphs. It is related to the kernel K' defined in the same way with 
the dressed propagators 

G = G'° + G'^K'G (7) 
where G'" is the product of the dressed propagators 

Gf = ^- for a fermion, (8) 

7j • Pi - TTii - + le 

Gf = —TT- — for a boson, (9) 

and Si the renormalized self-energy function. The transfer of the self-energies 
from the propagator to the kernel is achieved with [ pO| , p2[ 

K = K' + S = (G")-i - {G'^)-\ (10) 

In the two-body problem K' is simply the sum of all irreducible Feynman 
graphs. In the three-body problem, it is 

K' = {G'^r' + K',, (Gf + K'^, (Gi°)-i + K[,^ , (11) 

where K[j is the sum of the two-body (ij) irreducible Feynman graphs while 
-^123 is tfi^ sum of the three-body connected irreducible Feynman graphs. In 
the four-body problem and beyond, there appear mutually unconnected kernels 
in the expression of K' (like K'12 and K!^^ for example) , and we shall have to 
be careful to avoid overcountings in the expansion of G (see next section). 
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In ref. |Q, we built a 3D reduction of the two-fermion Bethe-Salpeter equa- 
tion around an instantaneous approximation of the Bethe-Salpeter kernel, the 
difference with the exact kernel being corrected in a series of higher-order contri- 
butions to the final 3D potential. This potential being not explicitly hermitian, 
we performed a second symmetrizing series expansion at the 3D level. The result 
of the combination of both series turned out to be independent of the starting 
instantaneous approximation, and was also obtainable directly by starting with 
an approximation of the free propagator, based on integrals in the relative en- 
ergies instead of the more usual 5— constraints. Furthermore, the method was 
easily generalizable to a system of N particles, consisting of any combination 
of bosons and fermions. In the present work, we shall however present our 3D 
reduction in yet another way. Let us first define the 3D free propagator 

N 

dpoG^ipo) EE 6{Po-^pio)dpio...dpNQG°{piQ,...pNo) 

i=i 



(-2i7r) 



with 



Po-S + iePo 



rg'P (12) 



S = EiA+-A-), (13) 



A±=Af...A±, Af = ^^, r = A+ + (-)/+iA" , (14) 



E = Y,E^, Lu = 2''Ef+i...EM (15) 



for f ^ 0, so that 



rg° = — + (-)^+i . (16) 

^ Po-E + ie ' Po + E^ie 

When / = 0, (bosons only), we have no r and no P (we can replace them by 
1 in (^) and 

_ 1 1 _ 2E 

^ ^Po-E + ie Po + E -ie ^ Pg - E^ + ie' ^ ' 

Let us now define a 3D full propagator in the subspace, with rt?" as free 
limit: 



(-2i7r) 



\jv-i '^V^ J dp'ndpoG{p'„,po) f3Ty/uj. (18) 
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The integrations with respect to the relative energies preserve the positions of 
the bound state poles of G (it can also be shown that the physical particle- 
particle, particle-bound state or bound state-bound state scattering amplitudes 
are also preserved If we write G as 

G = G° + G°rG°, T = K{l-G^Ky^ (19) 

we get 

5 = rgO + Tff <T>Tg° (20) 

with the definition 

<A>= / / dp',dpoG°ip',)A{p'„po)G'{po)P^. (21) 

To the 3D full propagator g we can associate the 3D equation 

= gOrV^0 (22) 

where V is such that 

5-' = r(g°)-i - K (23) 



The comparison with (20) gives 

<T> = V(l-rg"F)-i (24) 

or, conversely 

V =<T> {l+Tg'^<T>)-^ =< K{1-G°K)-^ > {l+Tg° < K{l-G°K)-^ >)-^ 

= < K{l-G°K)-\l+ > rg°< X(l-G°i^)-i)-i > 
= < K{1~G"K+ > Tg°<K)-^ >=< K{l-G^K)-^ > = < > (25) 
with 

= K{1 - G^K)-^ = K + KG^K + ■■■ (26) 
G^ G" - G^ G^ = > Tg° < (27) 
or, making the dependence on the relative energies explicit: 

G«(p^„Po) = G"(po)Mpo -Po) " G'{p'„po) (28) 

with 

G\p'^,Po) = g°(Po)/^ (_2j)i^-igO g°(Po)- (29) 

By substracting G^ from G'', we remove the leading terms which come from 
the residues of the positive-energy poles of G^ . This ensures in principle the 
decreasing of the terms of the series (p9). 
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The choice of g as the integral of G between two r operators defines 
our 3D reduction (in the two-body case it is also Phillips and Wallace's re- 
duction). A different approach, in the two-body case, the constraint approach, 
consists in defining < T > by fixing the relative-energy arguments of T{pq,pq), 
for example by putting one initial particle and one final particle on their positive- 
energy mass shell jl^. There is no problem with this approach for N = 2. For 
TV = 3, we should put two initial and two final particles on their positive-energy 
mass shell. For the two-body terms, this would give a constraint too much. This 
difficulty can be avoided by applying different constraints to the different terms 
of T |2^. This leads to a set of coupled equations which can not, without 
approximations, be reduced to a single equation. 

In we derived equation (^) from the inhomogeneous Bethe-Salpeter 
equation (H), with the definition 



(30) 



A simpler reduction method, which we shall adopt from now, can be obtained 
by keeping only the positive-energy part A+ of r. In this case, the equation 
obtained after the reduction of the Dirac spinors into Pauli spinors will be pOt] 



f 

n 



2E,, 



f 

n 

i=l 



2E, 



E,, 



if. 



(31) 



where v is the large-large part of V. This simpler reduction method can be 
defined by the choice of g as the integral of G, now between two positive- 
energy projectors A"*" (instead of r in (|8|)). This 3D full propagator is, in 
configuration space, the retarded part of G taken at equal times. It has been 
the starting point of a 3D reduction by Logunov and Tavkhelidze Q (in the 
two-fermion case), followed by Kvinikhidze and Stoyanov pTf (three fermions) 
and Khvedclidze and Kvinikhidze ||2^ (four fermions). 

3 The four-body problem. 

In this section 3, we shall neglect, for simplicity, the radiative corrections (omit- 
ting thus the E's and all the primes in the expression of A'). 



3.1 Bethe-Salpeter equation for four particles. 

For four particles, we have ^ p7| 

Go y^O /^O /^O 

— Lt-|^ Lzg Lz^ 

K = A'i2.34 + Ari3.24 + Ari4.23 



(32) 
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+ ifl234, (33) 

with 

i^l2,34 = Ki2 (G§ G° + if 34 (G° G° - if 12 if 34 ( 34) 

and similarly for if 13, 24 and if 14, 23- The three lines of (p3) contain the two, 
three and four-body kernels respectively. The last term of ( |34| ) is a counterterm 
which has no equivalent in the two and three-body problems. If we write the 
expansion of G 

G = G" G°if G° G°if G"if G" + • • • (35) 
and collect the terms containing one if 12 with one if 34, we get indeed 
G = ■ ■ ■ - G"Ki2K3iG" + G"ifi2(G!]G*l)"^G"if34(G?G^)"^G° 

+ G°if34(G; G° )-iG°ifi2(G§ G° r^G" + ■■■ (36) 

Since if 12 and if34 commute mutually, the two last terms of (^6|) correspond 
to the same graph which appears thus twice, while the first term is again the 
same graph with a minus sign. This mechanism works at all orders. The kernel 
-F'^12,34 is indeed given by 

i^l2,34 = (G° - (Gi2,34)"' (37) 

where Gi2,34 is the full propagator obtained by keeping only the interactions 
if 12 and if 34. It factorizes into 

(Gl2,34)"^ = (Gi2)-1(G34)-^ = [(G°G°)-l-ifi2] [ (G^] G^ ) " ^ - if 34 ] 

(38) 

so that ( |37| ) gives (^). This problem of commutating kernels was not present 
in the two and three-body problems. It has been easily overcomed here in the 
four-body problem. It will become much more complicated in the five-body 
problem and beyond |2^, . 

3.2 3D reduction of the four-body Bethe-Salpeter equa- 
tion. 

For definiteness we shall work on the Bethe-Salpeter equation for four fermions 
and perform the reduction based on the positive-energy part of g^. Then: 

A+ 5" = — (39) 

<A> = ^_A+(Po-i?) / dp'odpoG°{p'o)A{p',,po)G°{po)f3A+{Po-E) 



(40) 



8 



and the 3D equation will be 

(/. = gV</) (41) 

with 

V=<K'^> = <K> + < KG'^K >+■■■ 

^<K> + < K{G°- > g" <)K > + • • • 

= <K> + < KG^K >-<K>g"<K>+--- (42) 

We shall now compute the two and four- vertex terms, keeping only the two- 
fermion kernels: 

< X'J^ >(4) ^ < /s:i2,34 > + • • • (3 terms) 

+ [ < Ki2,34 G"Ki4,23 > - < if 12,34 > g" < ^^14,23 >]+••• (6 terms) 
+ [ < Ki2M G° A'i2,34 > - < is:i2,34 > 5° < 12,34 > ] + ' ' ' (3 terms) (43) 
where 

< ifl2,34 >-< {Gl Gl )-!> + < if34 (G? Gl )-^>-< Ki2 K34 > 

(44) 

< ifi2,34 G°Xi4,23 >^"^ - < K,2 {G^ G^ )-'G°K23 (G? G^ )'' > 

+ < K^2 [Gl Gl )-'G'>Kii (G° G" )-'>+< X34 (G? G" )-ig"A'23 (G? G° )-' > 
+ < X34 (G? G° )-iG°Xi4 (G° G° )-i > (45) 

< ifi2,34 G0Xi2,34 >("^ = < K,2 (G° G° )-iG0ifi2 (Gg G° )-i > 

+ < (G? G° )-iG0if34 (G; G° )-!>+< A-12 (G° G^; )-'G°K34 (G? G^ )-i > 

+ < X34 (G? G^ )-'G°Ki2 (GO G^ )-i > (46) 

The various terms of (|4^)-(^6|) are represented in figure 1. The two last terms 
of (|4^ ) are equal to < K12 K34 > . In the expansion of G this term will thus 
appear only once, with a plus sign. 

Let us now examine the counterterms ~<K>g'^<K>, which could 
be represented by the 8 last graphs of figure 1, with a vertical line separat- 
ing them into two parts. In the two and three-particle problems we have a 
cancellation between the corresponding contributions of < KG^K > and of 
— <K>g^<K>, when the two-body kernels are approximated by (rel- 
ative and total) energy independent kernels and the propagators by positive- 
energy propagators. This cancellation of the leading terms leads in principle 
to a decreasing of the contributions of the higher-order terms of the series 

< K > + < KG^K > H . Wc shall check this cancellation here in ( p| ) 

by assuming that the two-body kernels contain positive-energy projectors and 
are independent on the energies ( Kij = A'j Kij Pif3jA'^ A'^' and is independent 
of paij, Poij ~ we shall speak of the "positive-energy instantaneous approxima- 
tion" , although " instantaneous" usually refers only to the independence on the 
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two-body relative energy) . It is then easy to verify that the cancellation occurs 
for the 4th to the 9th graphs of figure 1. We remain thus with 

< >(4) ^ < (GO Gl )-!> + < X34 (G? > + • • • 



< K3i> - < Ki2 (G° Gl > 5° < Ksi (G? G! 

^0 r^O \-l 



■0 /-<0 -j-l 



> 

< X34 (G? G° > 5° < i^i2 (G§ G° > + • • • (47) 



The two-vertex contributions become simply the sum of the six two-fermion 
potentials: 

V = V12 + V3i + Vi3 + V24 + Vii + V23 , 



(48) 
(49) 



For < K12 K34^ >, we have 
< K12 K34 > = 

(Po - E') 



-1 
2i7r 



1 



dPi2Q dP3A0 5 {Po ~ P120 

1 



P: 



340 ; 



Pl2 



E'2 



P3An — En — E, 



Pr. 
1 



34 



-f340 ^ E3 



E4 + ie 



El- E2 + ie 
(Po - E) 



(50) 



^3 -^4 ^ 

Let us perform the integration with respect to P120 and close the integration 
path clockwise. We have to consider the poles at P120 = i?i+i?2 and P120 = 
E'l + E'^. We obtain 



< K12 K34 > : 
+ V12 



Pn-E' 



El + E2 — E[ — E'2 
1 



Vi: 



1 



^^34 



Po — El — E2 — E!^ — i?4 



34 



E[+E'2-Ei-E2 ■ Po - S( - - S3 - ^^4 ■ 
Writing then 

Po-E' = (Po - Si - ^2 - S3 - ^4 ) + (El +E2-E[-E'2) 
Po - E = [Po- E'l- E'2- E3- E4) + [E'l +E'2-Ei-E2 



(51) 

(52) 
(53) 



in ( pT| ) , we see that the contributions of the first terms cancel mutually, so that 
we remain with 



< K12K34 V12V3A 



1 



Pq — El — E2 — E'^ ~ E'4 Pf) — E'^ — E2 — E3 — E4 

(54) 

These two terms are cancelled by the two last terms of ( ^ ) respectively. 
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We have thus seen that, in the case of instantaneous positive-energy kernels, 
the contribution of the four- vertex terms vanishes, while the contribution of the 
two-vertex terms is simply the sum of the six two-fermion potentials. Let us 
now examine a typical set of six- vertex graphs, the ones which contain two K12 
with one K^^. Their contributions to the potential are represented in figure 
2, with a vertical line when G" is to be replaced by > < . After the mu- 
tual cancellations of some identical contributions, we obtain the corresponding 
contribution to < T >, minus the four graphs containing one > 5*^ <, plus 
the three graphs containing two > < . Let us now consider again the case 
of instantaneous potentials with positive-energy propagators. We have shown 
above that, in a sequence K12 K^^, we must consider the two possible orders. 
Here, in a sequence K12 K12 K^i, it can be shown that we must consider the 
three possible orders, so that the first graph of figure 2, e.g., becomes equal to 
the sum of the three last ones. In figure 3, we expand the graphs of figure 2 in 
the case of instantaneous kernels with positive-energy propagators, and we see 
that the sum of the resulting graphs is zero. 

4 Bypassing the Bethe-Salpeter equation. 

In section 3, we presented a 3D reduction method for the N-body Bethe-Salpeter 
equation. The starting homogeneous equation was written in terms of a kernel 
which was designed to reproduce the full propagator G by iterations of the 
inhomogeneous equation. We have seen that the writing of this kernel K was 
straightforward for N = 2, 3, less straightforward for = 4, and increasingly 
complicated for iV > 5. Our investigations on the iV = 4 case in section 3 
suggests us a possible way of avoiding the explicit writing of the Bethe-Salpeter 
equation. We saw in section 2 that the 3D potential is given by 

V =<T{l + G^T)-^ > = <K{1-G"K)-^ > (55) 

where G^ = G'~' — G^, while K contain counter-terms in order to avoid the 
apparition of topologically identical diagrams in the expansion of T = K{1 — 
G^K)~^. It is however possible to adopt a simpler algorithm for T : 

T = [K'^'il - G"if '''' . (56) 

By K^^- we denote the sum of the irreducible interactions (as (|3^^4|) without 
the three counter-terms) and by SG (single counting), we mean that the di- 
agrams which appear two or more times in the expansion of T must be kept 
only once. The 3D potential will then be written 

V = < K'^{1 - G"K^^)-^ >^'^ (57) 

in which, after the expansion of the series and the splitting of G^ into G'^—G^, 
we shall remove the duplicating diagrams. This makes the writing of the 3D 
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potential as straightforward for iV>4 as it was for A^= 1 or 2. As an example 
in the four-particle case: 

(Xi2G^K34 + K^AG''Ki2f'' 

= Ki2G°K34 - Ki^G'Ksi - K34G'Ki2. (58) 

With positive-energy instantaneous interactions, the term in G^ will be can- 
celled by the two terms in G^ . 

The writing of the Bethe-Salpeter kernel could be done in a quite similar 
way: 

K = {T(l+ G^T)-i}^,^^„ = {[i^^«(l-[G°-G^]i^^«)-i]"^}^^ 

^(59) 

Here, G^ is simply a temporary renaming of G*^, which indicates that terms 
like AG^ B and BG^ A must always be kept both. We know that K, unlike 

contain only a finite number of parts For iV = 4, e.g., eqs. (|3^p3) 

show that K contains 11 subkernels (groups of irreducible graphs) in 
plus 3 counterterms, coming from the second order in K^^ (we get one of these 
counterterms by replacing G^ by G° in (psf)). 

The operator T is directly given by the Feynman graphs if we neglect the 
radiative corrections. If not, we can use 

G° + G°TG" = G ^ G'° + G'"T'G'° (60) 

with 

G'° = G"(l-I]G")-i = (l-S^G^)-^ (61) 

i 

The operators G and T' are directly given by the Feynman graphs. For T, 
we have 

T = S(1-G°S)"^ + {l-T,G"y^T' {1-G"T,)-\ (62) 



5 Conclusions 

In a previous work we were in search of a 3D reduction method for the 
two-fermion Bethe-Salpeter equation based on an unspecified positive-energy in- 
stantaneous approximation of the Bethe-Salpeter kernel. We performed a series 
expansion around this approximation, followed by an integration on the relative 
energy and a second series expansion, at the 3D level, in order to render the 
resulting 3D potential symmetric. After combining both series, we found that 
we had in fact built a kind of propagator-based reduction, using an integration 
with respect to the relative energy, in full contrast with the usual constraining 
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propagator-based reductions, which use constraints. Furthermore, this method 
was easily generalisable to systems consisting in any number of fermions and/or 

bosons. 

In the present work, the increasing difficulty of writing the Bethe-Salpeter 
kernel for N > 4 suggested us a direct way of writing the 3D potential without 
the need of first writing the Bethe-Salpeter kernel explicitly. This writing of the 
3D potential is straightforward and valid for all N, with the simple prescription 
of removing the duplicating graphs which appear when N > 4., 
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Figure 1. Typical terms of <K+KG K>. 
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Figure 3. Two (12)-one (34) terms of <KG KG K> 
in tiie positive-energy instantaneous approximation. 
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Figure 2. Xwo (12)-one (34) terms of <K:G KG K:>. 



